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' We discuss fixed target and ISR inclusive photon production and attempt a comparison between 
theory and experiments. The dependence of the theoretical predictions on the structure functions, 
and on the renormalization and factorization scales is investigated. The main result of this study is 
that the data cannot be simultaneously fitted with a single set of scales and structure functions. On 
the other hand, there is no need for an additional intrinsic /c^ to force the agreement between QCD 
predictions and experiments, with the possible exception of one data set. Since the data cover almost 
overlapping kinematical ranges this raises the question of consistency among data sets. A comparative 
discussion of some possible sources of experimental uncertainties is sketched. 
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Despite many years of intense experimental |, |, ||, |, |, |, 0, |, |, and theoretical [|ll], |l|, 0] 
efforts the inclusive production of prompt photons in hadronic collisions does not appear to be fully 
understood. No consensus has been reached concerning the phenomenology of these processes. An 
attractive possibility made use of the intrinsic ambiguities of fixed order perturbation theory to define 
the various unphysical scales entering the theoretical predictions by means of various "optimisation" 
prescriptions |15, 0, 17, |l^]: an excellent agreement [0, |l^, 20, 21] between theory and experiments 



over the whole available range of energy and transverse momentum was thus obtained with a single set 
of structure functions and a unique value of A^^^ . More recently however it has been proposed to fix 
the unphysical scales at some arbitrary "physical" values and to introduce an extra non-perturbative 
parameter, called the "intrinsic transverse momentum" of the partons in the hadrons which is fitted to 
the data at each energy [^, ^ |25[ . This parameter increases with energy (technically speaking 
it is therefore not an "intrinsic" momentum) and this is interpreted as taking into account the effects 
of multiple soft gluon emission associated to the hard partonic scattering. 

Motivated by the recent publication of two new sets of fixed target data and using the latest up- 
to-date theoretical calculations we discuss the phenomenology of inclusive prompt photon production. 
The new data are those of the UA6 proton-proton and antiproton-proton experiment at a center-of- 
mass energy ^/s = 24.3 GeV Q and of the E706 proton-Beryllium and pion-Beryllium experiment at 

= 31.6 GeV and at ^/s = 38.8 GeV To avoid further ambiguities associated to the criteria of 
isolation we do not discuss the colliders data on the production of isolated prompt photons. 

We first set the theoretical framework by recalling the main features of the complete next-to- 
leading order (NLO) calculations in perturbative Quantum Chromodynamics (QCD). The intrinsic 
uncertainties of the NLO expressions are related to the choice of three arbitrary scales: the renor- 
malization scale, the factorization scale associated to the initial state collinear singularities and the 
fragmentation scale related to the the collinear fragmentation of a parton into a photon. A rather 
complete numerical study of these uncertainties is carried out. The main feature is that there is no 
stability zone in the fixed target energy range unlike what has been observed at higher energies [ p!^ . 
Arbitrarily choosing the scales we show that we can get a reasonable agreement with all considered 
data except with the E706 data obtained using nuclear targets and the older R806 data from 
the ISR. This may raise the possibility of inconsistency among various data sets. 



1 Theoretical framework and ambiguities 

Because of the well-known anomalous photon component the cross section for inclusive photon pro- 
duction takes a more complicated form than for pure hadronic reactions. The differential cross section 
in transverse momentum Pj, and rapidity rj can be written as a sum of two components: 
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dpj,dr] dp^dr] dpj,dr] 

where we have distinguished the "direct" component a^" from the "bremsstrahlung" one a^'^°™. Each 
of these terms is known in the next-to-leading logarithmic approximation in QCD i.e. we have 
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where the parton densities in the initial hadrons -^i/hi find i^j/zij depend on the factorization scale 
M while the parton to photon fragmentation functions D^y/^ depend on the fragmentation scale M^. 
The renormalization scale ^ appears in the strong coupling as- The higher order correction terms 



to the direct and bremsstrahlung cross sections, Kf" ||T^, 14 1 and K^pj^ ||2g] respectively, are known 



and we shall use their expressions in the MS convention. The dependence of these functions on 
the kinematical variables xi, X2, z, ^/s,pJ. and r] has not been explicitly displayed. All the quantities 
entering the above equations have been either calculated {Kf" and -f^]'/™^) or have been determined 



(see e.g. [|19|, ^ for and ||3^, |31|] for D^/f.) at the required level of accuracy by next-to- 
leading order fits to the data. The knowledge of A , from deep-inelastic experiments for example, 

MS 

completely specifies the NLO expression of the running coupling OsifJ-)- It is important to stress that 
the scales //, M and Mp are arbitrary. The dependence in the first two scales partially compensates 
within each of eq. (^) and eq. (|^) between the lowest order term and the correction term. The 
dependence in the fragmentation scale is more complex because of the dual nature of the photon 
which acts either as a parton (anomalous part) or a hadron. We come back to this point in great 
detail below. For our present purposes, it is enough to know that the dependence of the photon 
fragmentation function is compensated partly by terms in K'^" of eq. (|2|) and partly by terms in if 
of eq. dH) . Since the scale is arbitrary, it is clear that eq. (§) and eq. dH) which depend essentially 
monotonically on the scale do not have an independent physical meaning: only the sum, eq. (Q), 
has a chance to be (relatively) independent of and therefore is interpreted as a physical quantity. 

The fundamental ambiguity of the perturbative prediction is due to the fact that the three scales 
/X, M and are not determined by the theory. Roughly speaking they are parameters which control 
how much of the higher order effects are resummed in as(/x), and -D^/fc respectively and how 

much is treated perturbatively in K^J'-' and K^p™. In order to illustrate this problem keeping however 
the numerical work at a reasonable level we proceed as follows. We fix arbitrarily the fragmentation 
scale Mp and the factorization scale M. The remaining scale /i is then determined by solving the 



algebraic equation |16, 18 1 



djda/dppd'n) ^ ^ 
dln{fj.) 

The result is illustrated in the left picture of fig. || for the inclusive production of photons in proton- 
beryllium reaction at ^/s = 38.8 GeV, p^, = 6.12 GeV/c. We see that at fixed M^ there is a stability 
point (minimum) of the cross section as M varies which is a feature of a reliable perturbative prediction. 
However this stability point does not go through an extremum when the fragmentation scale varies. 
An alternative choice for the renormalization scale fj, consists in setting, for example, fi = M rather 
than using eq. (Q). This leads to the right hand picture in fig. ^ again no stability is achieved 
in the fragmentation scale while a stable point (maximum) appears when varying fi = M at fixed 
> 24 GeV^. Therefore, we do not have a stable point on the M, M^ sheet and one is somewhat at 
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Figure 1: Variation of the inclusive photon production cross section for proton-beryllium at Ej^^^^ = 800 
GeV, Pj, = 6.12 GeV/c, r/ = 0., as a function of the factorization scale and the fragmentation scale; 
left; the renormalization scale is "optimized" according to eq. (^j; right; the renormalization scale is 
set equal to the factorization scale. 



a loss about which scale one should choose to make predictions. This pessimism should be tempered 
however by noticing that there is a wide domain of relative stability for the values < 10 GeV^ and 
< 300 GeV^ (left picture) or > 24 GeV^ where a maximum in M can be obtained (right hand 
picture). Despite the completely different scale variation patterns of the two pictures in fig. |^ it should 
be noted that, for a given fragmentation scale, the theoretical cross section estimates at the stable 
points are numerically very close: they differ by less than 1%. Finally, when increasing Ml from 24 
GeV^ to 300 GeV^ the optimal predictions decrease by less than 20%. It should be reminded that at 
much higher energies, in the Tevatron energy range and above, a stable point could be obtained for 



the inclusive photon cross section |13| 



To understand the reason of the present unstability under variation of the fragmentation scale one 
can make the following observations. Consider the evolution equation of the photon fragmentation 
function. It is given by (dropping the longitudinal variable z and denoting instead the convolution in 
z by the notation 0) 

dln(M2) - + 27r (5) 

^ 3=q,9 

where the P^j and Pji are the appropriate Altarelli-Parisi splitting functions which admit a perturba- 
tive expansion in a^. For consistency we keep the first two terms in our analysis, namely 

p _ p(0) , as{M^) (1) 
and similarly for Pji. The solution to eq. (^) can be written as the sum of the anomalous part and 
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the non-perturbative part: 

, , , , „NP 



D^/k{z,M^)=D^^^{z,M^) + D^/,,{z,M^). (7) 



AN 

The component D^j^i^z, Mp) is fully calculable and it is a particular solution of the full inhomogeneous 
set of eq. (|5|). The non-perturbative part (modeled following the Vector Meson Dominance (VMD) 
ideas) D^^j^{z, M^) obeys the usual homogeneous evolution equations of hadronic structure functions 

i.e. eq. (^) with P^j = 0. The Mp scale dependence of D^J'^^, appearing in du'^''''™, compensates 
between the two terms of the right hand side of eq. (P) as it should be for purely hadronic reactions. 

AN 

In contrast, the scale variation of -C^/^ associated to the inhomogeneous term P^^ is compensated 
by a similar term in K!^" of eq. (|2|). It turns out that, in the kinematical region probed by fixed 
target experiments, such a scale compensation is not properly achieved. Indeed, due to steeply falling 
partonic cross sections, the average value of the variable z in the bremsstrahlung process is rather 
high, z .8. In that domain, we can safely neglect the gluon fragmentation function (as well as the 
non-perturbative VMD input) and check that the approximation 

D,,,{z,Mp) = ^2^ + d^^(z) 



as{Mp 

b c^,{z) \n{Ml/K'_) + d.„{z), (8) 



F J 



with the moments (c(n) = /q^ dzz" ^c{^z)) of the c-yg(z) function given by 



= , JoV; ' (9) 



satisfy eq. @ (neglecting 0{as) terms). In the above equations the parameter b is the well known 
positive constant governing the evolution of the strong coupling in the leading logarithmic approxi- 
mation 



dln(M2) 



as{Mp)Y {l + bias{Mp)). (10) 



At large z, the function b c^q{z) is much smaller than the function 'P^q{z) which appears as the 

coefficient of the In(M^) dependence in eq. (§). This can be guessed from eq. (^) since Pqq^ {ri^ is 
negative at large n, equivalently large z, as required by the negative scaling violations of the quark 
fragmentation function. The relative smallness of b c^q{z) is understood by the fact that it represents 
the shape of the fragmentation function of the quark into a photon after multiple gluon emission 
while P^q (z) ~ (1 -|- (1 — z)'^)/z describes the fragmentation at the lowest order level (without gluon 
emission). The overall scale Mp dependence of eq. (|l]) is then a weighted average over large z values 
of the combination 

{bc,q{z)-P(^^{z))ln{Mp) (11) 

which is always negative in the large z range, in agreement with the results of fig. |l[ This point is 
illustrated more precisely on fig. ^ which compares, at a fixed value of Mp, the variation in z of D^/g 
of eq. |8| with that of the compensating term proportional to P-yq- As anticipated, the compensating 
term has a less steep dependence and largely dominates at large z. Since the effective value of z in 
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Figure 2: solid line; z dependence of the perturbative component of the photon fragmentation func- 
tion eg. dashed line; z dependence of the compensating term proportional to P^q (z)). The scale 
isM^=6 GeV. 



inclusive photon production at E706 energies increases from z ~ .85 to about z ~ .93 when p^, varies 
from 5 GeV/c to 10 GeV/c it is clear that expression ( |TT1 ) is negative in the region of interest^. 

This discussion concentrated only on the scale variation driven by the inhomogeneous term (P-yq) 
in eq. (^). The homogeneous term, included in our calculation, somewhat weakens the dependence in 
expression ( pT| ) but does not change the overall behaviour. We do not see how to overcome this basic 
instability of the perturbative prediction for the inclusive production of direct photon in hadronic col- 
lisions. Notice that the perturbative approach is not spoiled by the presence of large ln(l — z) terms, 
for they cancel each other (at order a^) between the direct contribution eq. (|2|) and the bremsstrahlung 
contribution eq. (0) 



yj. Recently, two groups [32| have studied, in hadronic processes, the resumma- 
tion of collinear and soft logarithms which become large in the limit where xi and X2 are close to 1. 
Although the average values of xi and X2 in eqs. (|2|,^) are of the order of Xj, , the resummed expressions 
might slightly change the values of the predictions obtained in this paper. The stability with respect 
to the scales /i and M could also be improved. These types of corrections are not considered here. It 
would be interesting to investigate their impact on the prompt photon phenomenology. 

We propose to bracket the theoretical results by chosing two extreme values of the fragmentation 
scale. Let us note that making the choice Mp = fi may accidentally reduce the overall scale dependence 
of the cross section since, for example, one power of as(/i) in eq. @ is compensated by the l/as(/i) 
factor of Dj/q- We do not advocate this choice although we will often use it for practical reasons in 
the following. 

Before closing these theoretical preliminaries, let us make a technical comment. It concerns the 
definition of as(/u). The solution of the NLO evolution equation of the coupling (eq. (p!0|)) leads to the 



following relation in the MS scheme 

b ln(/iVA^) = ^- + 6iln 

MS Ctsl/J-j 

^As increases the relative importance of the fragmentation terms decreases compared to direct photon production 
terms. 



1 + biasin) 



(12) 
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Collaboration 


Reaction 


Vs 

[(jG V J 


Pj, range 

V / CJ 


/rapidity 
range 


Xrp range 


WA70 




P P 


23.0 


p^ > 4.0 


-.35 < Xp < .45 


Xj, ^ .35 


UAeg] 


P P 


24.3 


Pj, > 4.1 


-0.1 < ?7 < 0.9 


Xj, > .34 


UA6|] 


PP 


24.3 


Pj, > 4.1 


-0.1 < 7? < 0.9 


Xj, > .34 


E706[|] 


p Be 


31.6 


Pj, > 3.5 


-.75 < ?7 < .75 


> .22 


E706[|] 


p Be 


38.8 


Pj, > 3.5 


— 1.0 < 7? < 0.5 


Xy > .18 


R806[ 


1 


P P 


63. 


Pj, > 3.5 


-0.2 < 7? < 0.2 


Xy > .11 


RllO]^] 


p p 


63. 


Pj, > 4.5 


-0.8 < 7? < 0.8 


Xj, > .14 


AFS/R807[| 


PP 


63. 


p^ > 4.5 


-0.7 < 77 < 0.7 


Xy > .14 



Table 1: Summary of the main features of the data sets. The transverse and longitudinal variables 
defined by Xj, = 2pj,/^/s and Xp = 2p^l \fs. 

where 6 = (33 - 2Af/)/127r and 61 = (153 - 192iV/)/(27r(33 - 2iV/)) (with = 4, the number of 
flavors). To obtain as{fj.) explicitly one may invert exactly (numerically) this equation, or one may 
invert it approximately (analytically) keeping the first two terms in the expansion of as(^) in terms 
of l/ln(;U^/A^ ). In the region of interest for this study the difference between the obtained values 

MS 

of as{fi) for a given fi may reach 5%. We advocate, and we will use in the following, the "exact" 
estimate of as(/i). 

2 Phenomenology of prompt photon production 

Let us first summarize the data sets used in the following discussion. This is done in table |^. To 
avoid uncertainties due to the pion structure functions we restrict our discussion to the case of proton 
induced reactions. Concerning E706 and the possible nuclear effects on the parton distributions in 
the Beryllium target we follow the same prescription as the collaboration and adjust the theoretical 
predictions with the factor A°'~^ with a = 1.04. In this way the resulting cross section is normalized 
per nucleon. The inclusive differential cross sections are shown in fig. || where they are plotted as a 
function oipj,. Already, we anticipate some difficulties when comparing the data sets: for example, the 
relatively large difference between the WA70 and UA6 cross sections at low p^ cannot be accounted for 
by the small ^/s change between experiments nor by the different rapidity coverage; one also observes 
almost an order of magnitude difference between the \JA6{pp) data and the E706(530 GeV) for an 
energy change from = 24 GeV to = 31.6 GeV while the E706(800 GeV) and R806 data at low 
Pp have almost the same size despite an energy increase from ^/s = 38.8 GeV to ^/s = 63. GeV. 

In the numerical calculations we will use several sets of recent next-to-leading order structure 
functions in the MS convention. CTEQ p7| ] and MRS |^] will be our main sets. However, in order 
to compare the present results with the "old" fits done using the WA70 data we will also make use of 
the original ABFOW |jl9| sets |^ Concerning the higher order calculations they have been performed 
using four light flavors and therefore 7?T-chaim = in the partonic cross sections and the higher order 
terms. However the threshold effects due to the charm quark mass are taken into account in the 
fragmentation function D^/c in such a way that at a scale Mp ~ TTT-charm one recovers all the correct 
logarithmic factors of the exact massive calculation |3^. The most recent parametrization of the NLO 

^In our previous phenomenological studies ^] the higher order calculations to the bremsstrahlung cross section 
were not available and a stability in the fragmentation scale could not be expected. That scale was then chosen to be s. 
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Figure 3: The differential inclusive photon cross sections as a function of transverse momentum. 
The cross sections are averaged over the x ^ /rapidity ranges shown in table Data from the NA24 
collaboration ^ are averaged in the rapidity range —.65 < rj < .52. 



parton to photon fragmentation functions of ref. [31| will be used. In this work, two parametrizations of 
the fragmentation functions -D^/j, based on fits to e+e~ data, are given: we use set II parametrization 
which corresponds to the large non-perturbative gluon input. In fact, little difference is observed 
in inclusive photon production in hadronic collisions between set I and set II predictions since, as 
discussed in the previous section, the effective fragmentation variable z is large and the fragmentation 
functions are then dominated by their perturbative components. 

We first discuss the comparison between theory and experiments for proton-proton and proton- 
antiproton (when applicable) collisions, using a fixed common scale jj, = M = Mp. Then we shall 
relax this constraint and study optimised predictions in the /i and M scales for a range of fragmen- 
tation scale in order to propose an estimate of the errors on the perturbative predictions. In the 
theoretical predictions we always include the NLO corrections to both the direct component and to 
the bremsstrahlung component (denoted ACFGP in the figures). 



In order to normalize our present results to the previous ones appropriate to WA70 data we start 
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our analysis using the ABFOW parametrization (A = 230 MeV, x G{x, Ql = 2 GeV^) = (1 - x)^). 

MS 

We get with all scales in the problem set equal to a good approximation of the old optimised 

results and also a good agreement with the WA70 data (x^ = 8.2 for 7 points using statistical errors 
onlyP]). Comparing with the new UA6 measurements, there is an excellent agreement with the data 
(see fig. ^) on the difference pp — pp (which is insensitive to the gluon density in the proton and to the 
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Figure 4: Comparison of the UA6 pp — pp data ^ with the NLO theory for various sets of structure 
functions using the common scale fi = M = Mp = p^jl. Only statistical errors are shown. 



parton fragmentation functions into a photon) while the pp data tend to be higher than the theory 
by roughly 20% and the pp data by almost a factor 2 for Xj, < 0.4 (x^ = 64.4/9). The normalized 
ratios data /theory are shown in fig. ^. At large the WA70 and UA6 data are mutually compatible 
within the error bars while there may be some relative normalization problem for the points below 
Xj, = 0.4. Turning now to the ISR results one notices that they are located at small Xj, values and 
that, for the RllO and R807 data, the ratio data/theory is flat as a function of x^,; furthermore it 
is compatible with 1 {x^ = 6.8/7 for RllO and = 13.7/11 for R8O70). In contrast, the R806 data 
show a marked decrease from 1.5 to 0.5 as Xj, increases. It is clear that the problem of the R806 data 
is not related so much to the overall normalization as to the slope of the Xj, distribution. In fact, it 
was noted before that the Pj, dependence of the R806 data was incompatible with the slope of the 
other direct photon experiments. Besides, it was not possible to obtain a good fit simultaneously to 
the deep-inelastic data and to the R806 data [19|. However all ISR sets, up to Xj, = 0.3 are point to 
point compatible taking into account statistical as well as systematic error bars. 



^As in the lowest pr point is not included in the evaluation, due to important systematic errors . 
"AU the above values are obtained using statistical errors only except in the case of AFS/R807 where only the 
sum of statistical errors and systematic errors added in quadrature is available. 
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Figure 5: Comparison of data with theory, normalized to the theoretical results, for the various exper- 
iments as a function of x^. The ABFOW structure functions are used and all scales are set equal to 
Statistical error bars are shown as full lines while systematic uncertainties are added in quadra- 
ture and shown as dotted lines. An extra 5% uncertainty on the energy scale is not plotted for RllO; 
for AFS and E706 data, statistical and systematic errors are shown combined in quadrature only. 



The main surprise comes when comparing the E706 data to theory: at both energies we observe 
that the experimental points are 2 to 3 times larger than the theoretical predictions at high Xj, and 
the ratio experiment/theory keeps increasing up to a value of 5 as x^, decreases. This had already 
been stressed by the E706 collaboration Q. 

A similar study can be conducted with the common scale in the calculation set equal to Pj,/3 
instead oipj./2. The corresponding results are displayed in fig. ^ Good agreement is achieved for the 
difference pp — pp and the UA6 pp data (x^ = 14.7/9). Agreement between theory and experiment 
is much improved in the case the UA6 pp data (x^ = 23.1/8) while, on the contrary, the comparison 
between WA70, RllO and R807 is not so satisfactory (all values are roughly multiplied by a factor 
3 compared to the analysis with scales Pj,/2). In fact, all these data sets fall somewhat below the 
theoretical predictions. Concerning E706, essentially the same situation as before prevails since the 
theoretical predictions are rather stable (they increase by about 20 %) when changing the common 
scale from Pj,/2 to Pj,/3, namely the theory still underestimates the data by roughly a factor 2 at 
large and much more at small x^,. 



The provisional conclusion to be drawn is that there is room for a minor incompatibility in the 
normalization of the UA6 and WA70 pp data at low Pj. while they are in perfect agreement at large 
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Figure 6: Comparison of data with theory, normalized to the theoretical results, for the various exper- 
iments as a function of x^. The ABFOW structure functions are used and all scales are set equal to 
Data points at scales such that (pj,/3)^ < Q"^ = 2. GeV^ are not included, which explains why 
the points at lower Xj. values in some experiments are not displayed. See fig. |^ for further comments. 



transverse momentum, or to put it differently, there may be some incompatibility in the p^, dependence 
of the two sets of data. The ISR data are also compatible with WA70/UA6 while the E706 results 
stand much above all other experiments. It is to be noted that the UA6 pp data are in good agreement 
with the theoretical predictions both in shape and normalization. 

We repeat our numerical exercises using the CTEQ4A1, CTEQ4M and MRS-98 structure functions. 
In the case of CTEQ4A1 we obtain very much the same phenomenology as before. As illustrations, 
we show the comparison between data and theory using CTEQ4M with all scales equal to Pj,/2 in 
fig. 7^ and using MRS-98-2 with all scales equal to Pj,/S in fig. ^. Fig. |^ looks very similar to the 



case of ABFOW with the smaller common scale Pj,/3: this is understood because of the larger value 
of A_ = 296 MeV for CTEQ4M (instead of 230 MeV) which yields a larger value of Og for a given 
scale and therefore larger theoretical predictions. The same comment would apply to the results based 
on MRS-98-2 with the common scale Pj,/2. Prom fig. there appears a perfect agreement between 
theory and the UA6 pp and pp data (the values are 3.12/8 and 5.5/9 respectively) and a tendency 
to reduce the disagreement with E706 to an overall normalization factor of roughly 50 % at the cost 
of overestimating the ISR data and the low p^, WA70 data. Let us remark that the overall improved 
agreement between theory and experiment is, in part, due to dropping a few experimental points at 
low Pj, because of the scale limitation. 

**The difference on the ratios data/theory observed for the E706 points, at low p^, between fig. ^and fig. 14 of ref. Q 
arises partly because in the latter work the higher order terms in eq. were not taken into account. 



10 



6 



o 

CD 



a 
o 



4- — 



3 



□ 
▲ 

O 



2 



O 




WA70 pp 
UA6 pp 
E706 pBe, 
E706 pBe, 
UA6 pp 
R 1 1 O pp 
RS06 pp 
AFS pp 



'530 
'BOO 



NLL ACFGP 
M=/j,= MF=pt/2 
CTEQ4M A = 296 MeV 
frag BFG 



O. 1 



0.2 



0.3 



0.5 



C) A5 



O. 



0.8 



Figure 7: Dependence in of the experimental cross sections normalized to theoretical predictions 
using the CTEQ4M distributions. All scales in the calculations have been set to Pj./2 and only data 
points satisfying the condition {p^/2)'^ > Qq = 2.56 GeV^ are kept. See fig. ^ for further comments. 



It appears difficult to reconcile the E706 data with the other fixed target data: clearly their 
normahzation is higher. At large x^. values {Xj, > .3), one can roughly estimate from the previous 
figures that, when normalized to any perturbative NLO predictions, they are above the other data 
sets by a factor ranging from at least 50% (compared to UA6, all scales p^./'i) to 2.0 (compared to 
WA70, ah scales Pj,/2). 

A further striking fact, which is obvious on all figures, is the rise of the E706 cross sections compared 
to the theory as x^ decreases below 0.25. We point out here that this rise is not systematic as claimed 



in |22|. In fact, it is not observed at all on the RllO and R807 data and it is rather weak on the R806 
data which cover a lower x^ range than E706. Such a conclusion was also reached by the authors 
of 

In conclusion, no set of structure functions is able to accommodate the recent E706 data despite 
allowing a reasonable variation in the choice of scales. In contrast, the NLO calculations are in 
agreement (within the rather large experimental and theoretical error bars) with the other data sets 
which bracket the E706 data in energy and which, taken together, cover the large Xj, range of E706. 



^^In Ig^l (see also |2l|]) a better theoretical input than in 
bremsstrahlung component, was used. 



221, in particular a complete NLO treatment of the 
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Figure 8: Dependence in x^. of the experimental cross sections normalized to theoretical predictions 
using the MRS-98-2 distributions. All scales in the calculation have been set to Pj./'i and only data 
points satisfying (p^/S)^ > Qq = 2. GeV^ are shown. The condition p^ > 5. GeV/c (see sec. ^ 
translates into x^ > .32 for E706 (530 GeV) and x^ > .26 for E706 (800 GeV) and Xj, > .16 for ISR 
data. See fig. for further comments. 



3 Further phenomenological considerations 

In order to bring their data in agreement with theory the E706 coUaboration advocates the use of 
an extra parameter identified as a measure of the transverse momentum fluctuations of the cohiding 
partons ||3^0. In the present phenomenological implementation of this effect both the normahzation 
and the shape of the p^ spectrum are affected so that the E706 collaboration is able to reproduce their 
data with values in the range 1.2-1.3 GeV/c. For UA6 the choice = 0.7 GeV/c gives a good 



agreement with the data [24|. From these values one would expect kj, > 1.5 GeV/c to be appropriate 
for ISR data which would destroy the rough agreement between data and theory displayed in the 
figures. From this we conclude that, including k^, effects may help some data sets (E706) to agree 
with theoretical predictions but it simultaneously destroys the agreement with other data sets (WA70, 



ISR) with theory |34]. Globally the phenomenology of "low" energy photon production is not improved, 
on the contrary! More precisely, except for the E706 data there is no need for an extra parameter to 
obtain agreement between data and theory. 

Some more comments can be added concerning "/c^" smearing effects. There is no definite theoret- 
ical method to parameterize such effects and as a result different groups obtain rather different shifts 



tt 



In practice, more than one parameter is needed to model this effect l22 
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of the differential cross sections, specially at low values of p^., as it has been briefly discussed in [p^ . 
Besides, the fitted value of "fcy" depends on the estimates used for the perturbative cross sections. 
In the low region, the perturbative predictions become notoriously unstable because no optimum 
in the factorization and renormalization scales is achieved. Therefore one observes large variations in 
the perturbative estimates under changes of scales. We comment on this point next. 
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Figure 9: Dependence in Xj, of the experimental cross sections normalized to theoretical predictions 
using the CTEQ4-M distributions. All scales in the calculations have been set to p^/2> and small p^ 
points have been dropped as explained in the text. 



We have checked in numerical studies similar to those leading to fig. || that no stable point, at 
fixed fragmentation scale, can be obtained for E706 at 800 GeV and Pj, < 5. GeV/c, for example. The 
three-dimensional plots show a monotonic variation under changes of scales as would be obtained in a 
leading logarithmic calculation. We would like to interpret this fact as an indication that perturbation 
theory is not valid in that range since variation of the arbitrary scales leads to large variation in the 
predictions. For quantitative comparison between theory and experiments such points should not be 
included. This phenomenon occurs for the WA70/UA6 data for p^, values below 4.2 GeV/c. Let us 
note that at these larger values of transverse momentum the perturbative predictions are stable against 
moderate k^, effects {< kj, > < 0.4 GeV/c). If we accept this reasonable limitation the disagreement 
between data sets and the theoretical predictions resides essentially in an overall normalization factor 
and not in the shape of the Xj, distributions as seen in fig. ^ and fig. |9[ 

To gauge the agreement between theory and experiment more quantitatively and in order to roughly 
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[GeVJ 


(stat.+syst.) 


(stat.+syst.) 


(stat.+syst.) 


y2 

(stat.+syst.) 


WA70 m 


P P 


23. 


9.1/7 


21. /5 


5.7/7 


37.6/7 


UA6 g 


p p 


24.3 


13.7/9 


2.34/6 


19.0/9 


3.12/8 


UA6 g 


p p 


24.3 


5.2/10 


7.45/7 


6.44/10 


5.51/9 


R806g] 


p p 


63. 


55.4/11 


91.6/11 


91.6/11 


80.1/11 


R110§ 


p p 


63. 


5.7/6 


14.1/6 


3.9/6 


11.2/6 


AFS/R807i] 


p p 


63. 


27.8/10 


71.9/10 


18.2/10 


71.6/10 



Table 2: Table of y^ values normalized to the number of experimental points. The statistical errors 
and systematic errors added in quadrature are used in the evaluation of values. Data points at 
scales lower than Qq are not included, which explains why the number of considered experimental 
points is smaller for the scale Pj,/S than for the scale Pj,/2. Qq = 2.56 GeV^ for CTEQ4M and 
Ql = 2. GeV^ for MRS-98-2. Only the experimental points with Pj, > 4.2 GeV/c for WA70/UA6 and 
with Pj. > 5. GeV/c for ISR are kept. For RllO a 5% uncertainty on the energy scale is not included 
in the systematic errors. 



take into account the systematic uncertainties we proceed in the following way. First, we calculate for 
each experiment the value per number of experimental points using the statistical and systematic 
errors added in quadrature. The results are shown in table |2| except for the E706 data for which the 
systematic errors are not yet tabulated. Both sets of structure functions (CTEQ4M and MRS-98-2) 
lead essentially to the same conclusions: no agreement with R806 but satisfactory agreement with 
the other experiments. More precisely, the UA6 pp data are in excellent agreement with both sets 
of structure functions, independently of the scale choice which makes them a very good channel to 
extract the value of A__. . The UA6 pp data coincide very well with both sets of predictions based 
on the scale choice Pj,/2> while clearly the other experiments prefer the larger scale p^/2 or even a 
larger scale (alternatively, a lower A_ value would fit the WA70, RllO and R807 data very well as 
can be seen from fig. ^). Turning now to the E706 data the ratio data/theory is essentially flat in the 
Xj, region considered as perturbatively reliable (see figs. ^ and |^ but the data are in excess of the 
theoretical prediction by at least 50% both at 530 GeV and 800 GeV. To summarize, it could be said 
that all the data sets considered could be made consistent with each other (and with the theory) if 
one would allow for rather large relative overall normalization shifts. 

In table ^ we show the result of such a study and give the y^ values based on statistical errors but 
fitting an overall normalization factor for each experiment. Clearly, one gets for all experiments (except 
R806) very good values at the expense of sometimes unusually large factors A. An interesting case 
is E706: at 530 GeV the small y^ values could have been anticipated from the figures while at 800 
GeV the rather large values displayed in the table are clearly due to the data points at small x^, with 
small statistical errors (but large systematic uncertainties). Obviously our naive procedure cannot 
handle such effects correctly. It would be worthwhile, in this respect, to analyse all data including the 
point to point as well as the overall systematic errors properly. We can just conclude from table ^ 
that the price of obtaining consistency among all experimental results as well as a good agreement 
between theory and experiments is to allow normalization shifts in the data which are sometimes well 
in excess of the allowed ranges estimated from adding the different systematic errors in quadrature. It 
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I'stat "1 
normaliz. 


WA70 g 




23. 


7.55/7 
A =1.46 
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A =1.27 
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A =1.21 
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A =1.61 


UA6 |g 




24.3 


5.55/9 
A =0.79 


3.4/6 
A = 0.98 


5.3/9 
A =0.74 


3.5/8 
A = 0.93 


UA6 § 


p p 


24.3 


6.8/10 
A =0.94 


5.5/7 
A =1.11 


5.9/10 
A =0.91 
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A =1.06 


E706 g 


p Be 


31.6 


15.7/13 
A =0.52 
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A =0.64 


17.8/13 
A =0.50 


15.1/13 
A =0.64 


E706@ 


p Be 


38.8 


64/13 
A =0.50 


51.9/13 
A =0.61 


59.7/13 
A = .48 


48.7/13 
A =0.60 


Rsoeg] 
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63. 


329/11 
A =0.99 


278/11 
A =1.17 


348/11 
A =0.93 


281/11 
A =1.12 


R110[|] 


VP 


63. 


1.91/6 
A =1.26 


1.62/6 
A =1.50 


1.88/6 
A =1.19 


1.61/6 
A = 1.43 


AFS/R807|6| 


p p 


63. 


2.6/10 
A = 1.35 


3.2/10 
A = 1.58 


2.6/10 
A = 1.28 


2.7/10 
A = 1.58 



Table 3: Table of values normalized to the number of experimental points. The statistical errors are 
used in the evaluation of values hut, for each experiment, an overall normalization factor is allowed 
to vary. Only points with Pj, > 4.2 GeV/c for WA70/UA6 and with Pj, > 5. GeV/c for E706/ISR are 
kept. Strictly speaking, A should set to 1 for AFS/R807 since the quoted errors include the systematics. 

would appear that CTEQ4M, with the common scale in the theoretical predictions set equal to p^j'i^ 
leads to the best phenomenology in the sense that the normalization factors A are closest to 1, still 
keeping small values. However, because of the smaller number of experimental points kept in the 
analysis it is not possible to argue decisively in favor of this set rather than MRS-98-2. 

It was mentioned above that the shape of the E706 data at 800 GeV was not in agreement with 
the theory at small x^,. This raises the question of how changes of scales affect the shape of the 
predictions. This is studied in fig. ^ where the theoretical results for various choices are compared 
to the reference case with fi = M = Mp = Pj./2. Two extreme cases are considered, = 2 GeV^ 
and = 300 GeV^, the other scales being optimised in the sense of eq. (^. The effect of changing 
the scales consists essentially in an overall shift in normalization but the increase in the cross section 
associated to the small Mp scale is not enough to bring the theory and the E706 data in agreement. 
A relative enhancement is obtained at small Xj, with the (extreme) choice = 2 GeV^ but it is 
not sufficient to account for the shape of the E706 / 800 GeV data at low x j, . A smaller spread in the 
theoretical results would, of course, be obtained had we used Mp scales proportional to Pj, . 

A final comment concerning the reliability of the theoretical predictions is related to the uncertainty 
in the behaviour of the quark distributions at large x. Recently, it has been pointed out that nuclear 
effects in the extraction of the quark distributions from deep-inelastic data on Deuterium had not 
been fully taken into account Taking these effects into account leads to a much larger d/u ratio 
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Figure 10: Study of the sensitivity of the cross section to changes of scales for E706 at 800 GeV. All 
results are normalized to the cross section calculated using ji = M = Mp = 



at X > 0.5, increasing with x. As a consequence, the predictions for some prompt photon production 
rates could be affected. In fact, pp reactions should not be modified because the photon production 
rate is dominated by Compton-like scattering and therefore it involves the same combination of quark 
distributions as probed in deep-inelastic scattering on a proton, namely (4ti + (i), which is of course not 
sensitive to the effects discussed in [SC]. Such is not the case for nuclear targets. For pN scattering. 



where N is an isoscalar nucleus, a^^ ~ 5 [u + d) and an increase 5d in the d quark distribution 
manifests itself by an increase 6a^^ ~ 156d/4 where the constraint 6u ~ —6d/4: has been used. In [^] 
the effect of nuclear binding is parametrized as an increase of 5{d/u) ~ (.1 it .001)(1 + x)x of the d/u 
ratio. This translates into an increase of the theoretical predictions for pN cross sections 

which could be as large as 20% when the relevant < x > is large enough. This is a very rough discussion 
since, if nuclear effects have to be disentangled when extracting the d distribution in deep-inelastic 
experiments they should be included (back) when making predictions for prompt photon production 
on (different) nuclear targets! This point certainly deserves further detailed studies. 

The predictions for pp could also be modified due to the important contribution of the annihilation 
channel qq — > 7G specially at large p^, where the increase in the d/u ratio should lead to a slight 
decrease of the cross section (a quick estimate leads to Sa^'P / a^'P\ann ~ —5{d/u)/2 for a large enough 
Pt)- 

For lack of parton distributions incorporating this effect we do not pursue a quantitative study of 
this problem. 
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4 Experimental data sets 



In view of possible incompatibilities between various data sets, let us review some specific features of 
the experiments. Details on systematics uncertainties are not always available. For the latest E706 
data, details on energy scale uncertainty are available. For other issues, we comment on the early data 
analysis which may not apply for the large statistical sample. It would be interesting to have such 
information on the large statistical sample for a detailed comparison between data sets. 

4.1 Hydrogen or heavier nuclei 

In fixed target experiments, the data mainly have been taken at CERN with hydrogen (with the 
exception of NA3 data on isoscalar Carbon target) while the published direct photon cross sections 
with heavier nuclei come mainly from the E706 collaboration @, ^. A detailed study of the nuclear 
corrections is not available. 

4.2 Isolation criteria 

All the comparisons with theory have assumed fully inclusive photon production. In fact, the prompt 
photon samples are always obtained with some isolation criteria, at least to remove electrons and 
charged hadrons in the sample. These cuts depend on the granularity of the detectors which is an 
important ingredient of the shower reconstruction algorithms. The data presented are usually corrected 
for this inefficiency. 

These cuts were sizable for the pioneering ISR experiments. In the RllO sample, where the cuts 
were no charged tracks within 25 cm of the shower and no additional electromagnetic shower within 
20 cm of the candidate photon, no correction was applied and the data are presented as biased against 
bremsstrahlung |^] . A proper theoretical treatment of such cuts is not available. The R806 [Q] and AFS 

data, on the other hand, are corrected for such effects. The correction is 1.15 to 1.25 for the later 
but no details are given on the Pj, dependence of the correction. Part of the observed discrepancies 
between the ISR experiments at low could be attributed to this effect. 

For fixed target data, two kinds of cuts are applied, one on a charged track extrapolating to the 
shower of the photon candidate to remove electrons and charged hadron showers, and one on the 
longitudinal energy deposition to remove hadronic showers in the sample. For WA70 and UA6, the 
cut is no charged track extrapolating within 5 cm, respectively 1.5 cm, of the shower center. For 
E706, no details are given for the present data, but earlier data were obtained without charged track 
extrapolating to 1 cm of the shower center [^. For UA6 and E706 (530), due to the difference in energy, 
these cuts cover about the same phase space while for WA70 the cut is more drastic. Cuts on the 
longitudinal energy deposition are implemented differently by WA70 (at the shower tail) and E706 (at 
the shower front) and are corrected for. These two kinds of cuts have more effects on bremsstrahlung 
than direct photon production, i.e. more effects at low p^,. Small differences on cross sections may 
result from the different cuts and the different ways their efficiencies are estimated. They are surely 
bounded by the expected size of the bremsstrahlung contribution (which should not exceed 20%). 
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4.3 Backgrounds 

The experimental candidate photon samples are always contaminated by sizable backgrounds. With 
high granularity calorimeters such as those used by WA70, UA6 and E706, the backgrounds are 
mainly from vr*^. The amount of background depends on the calorimeter, for example granularity, and 
on the strategies adopted to reconstruct the showers. It obviously depends on p^, on rapidity, and is 
important at low Pj, in any case. The background corrections are performed with simulation of the 
known backgrounds in the detector, the simulation being checked on the reconstructed vr'^'s. 

Typical background fractions in the candidate photon sample are of order 50-60% at Pj,= 4.25 
GeV/c to 20-30% at p^=5.5 GeV/c (WA70, central rapidity); 70 to 80% at p^ of 4-5 GeV/c (E706 
early data). For UA6 pp data, the corresponding figure is about 40% and depends on while for pp 
it depends weakly on p^, and amounts to about 50% |3^. Typical uncertainties on the cross section 
quoted are 16% at Pj,=4 GeV down to 8% at Pj, of 8 GeV/c for E706 early data. For WA70, the 
uncertainties were in the 10-25% range with an additional 15% upper limit from simulation. UA6 
quotes specifically the various independent uncertainties on simulation in the 1 to 6% range. 

Small errors on background evaluation induces large errors on the extracted signals. It is therefore 
very important to understand the background precisely. This point deserves futher study as it may 
be helpful to understand discrepancies between experiments. 



4.4 Comparison of data sets 

Without theoretical prejudice, one can only compare data in a similar energy range. For the CERN 
fixed target data from the WA70 (beam of 280 GeV/c [|]) and the NA24 (beam of 300 GeV/c) 
experiments, "the agreement (at y=0) appears to be well within systematics and even statistical 
uncertainty" On the other hand, at ISR with ^/s = 62.3 GeV, "normalizations were about 

1.5 to 2 in relative doubt" ||3J 



Although most of the ISR prompt photon data have been taken at v^=63 GeV [Q, |5|, |^, data at 
lower -y/i | are reported as j/tt^ by the R806 collaboration (at ^/s = 31 GeV, 45 GeV and 63 GeV 

|4^). One can tentatively combine the data on "f/ir^ |41| at ^/s = 31 GeV with vr'^ production cross 
sections (using the super-retracted data from table 5 of Ref. [^] at = 30.6 GeV) |: the resulting 
7 cross sections turn out to be, within large experimental errors, compatible with the E706 results in 
the Pj, range 4.0 — 6.0 GeV/c. It is worth mentioning that this is just a rough cross-check as the R806 
collaboration did not publish a direct photon cross section at ^/s = 31 GeV ^. 

We have not included these data on the data/theory plots as they do not help to solve the observed 
discrepancies on these plots between the rough agreement with perturbative theory of the bulk of 
ISR data at low Xj, and the strong disagreement of the E706 data at higher x^, : within the deduced 
experimental errors, these ratios for the data at ^/s = 31 GeV are compatible with the E706 data at 
^/s = 31.6 GeV in the x^. range 0.26 — 0.39 but they are also compatible with 1. Unfortunately we 

*As for the published ISR data used in this present work, the discrepancy is in the slope in Pj, rather than in the 
overall normalization. 

^Data on 7/all neutral clusters have been reported at a/s = 44.8 GeV and 62.4 GeV by the RllO collaboration ]4o| ]. 
At fixed , the 7 / all neutral clusters decreases with i/s |Q while 7 /tt" shows little y/s dependence |^] . 

*These tt" cross sections are in fair agreement with measurements of the RllO collaborations at = 31 GeV [^ . 

^As a check of this rough procedure, combining the same way data on 7/^" at -^s — 63 GeV ^] with the super- 
retracted data at y/s — 62.8 GeV of Ref. [M one gets photon cross sections statistically compatible with Ref. H. 
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therefore cannot conclude on ISR/E706 comparisons. 



4.5 Summary on experimental data sets 



Other effects, such as luminosity, trigger efficiencies, absolute energy scale, may affect the 7 cross 
section. However they also affect the vr'^ cross section. It is more relevant to discuss them in this 
context [44]. 



5 Conclusions 



It appears that the present inclusive prompt photon data at fixed target and ISR energies are incom- 
patible: normalizing the theory on one set, the extrapolation of the theory to a slightly different energy 
completely misses the corresponding data. It does not appear very instructive to hide this problem by 
introducing an extra parameter fitted to the data at each energy: the agreement with some data sets 
is improved while the agreement with other sets is destroyed. Besides, the actual value of the extra 
parameter depends on the value chosen for the scales used in the perturbative predictions. Since, 
anyhow, the perturbative predictions become unstable at low p^, where the "intrinsic" momentum 
contributes most, the extracted value of this parameter suffers from a large error. It seems more rea- 
sonable to us to only use data at reasonably large values. In this region, the various sets appear to 
show less important discrepancies. In any case, an effort should be done by experimentalists to extend 
the discussion of sec. § to better understand the various sources of systematic uncertainties. Perhaps 
RHIC, with a center-of-mass energy as low as 50 GeV, could help resolve the important discrepancy 
observed between some data sets. 

As a final remark, let us notice that all direct photon experiments also measure inclusive vr'^ cross 
sections. It would be interesting to see whether the normalization discrepancies observed in direct 
photon experiments persist for the vr'^ cross sections. Here the theoretical predictions are less reliable: 
1) the cross sections depend on the fragmentation functions of quarks and gluons into tt"; 2) the 
effective value of the fragmentation variable z is in the range .8-. 9 and large ln(l — z) terms should be 
resummed; 3) the scale dependence of the cross section is large. For all these reasons, we do not expect 
a very good agreement between data and theory in the Pj,-range of fixed target experiments [^]. We 
can however use the theoretical cross sections as reference curves and investigate how data are situated 
with respect to the predictions. Preliminary studies show that fixed target data overshoot the 
theoretical curves by a few tens of percents, when CTEQ4M or MRST and the scales Pj./'i are used. 
Therefore all fixed target -k^ data display a similar behaviour compared to the theoretical cross sections 
and appear less scattered than the prompt photon data. This would indicate that the main problem 
rests with the extraction of the 7/7r'^ ratio. 
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